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Abstract 



We present the result for the first moment of the scalar and axial-vector 
current correlator in third order of the strong coupling constant ag and 
give the details of a recent evaluation of the pseudo-scalar correlator. The 
results can be used to reduce the theoretical uncertainty due to higher order 
corrections for the determination of fundamental parameters of QCD in the 
context of lattice calculations. 
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1 Introduction 



The computation of heavy quark current correlators in perturbation theory is of cen- 
tral interest for many phenomenological applications. For instance their low energy 
expansion can be related to moments, which have been determined with high preci- 
sion through the calculation of higher order corrections. Moments of the vector current 
correlator play an important role in combination with sum rules and the experimen- 
tally measured ratio R{s) = a{e'^e~ — >• hadrons)/cr(e"'"e~ — >• They can be 
used to perform a precise charm- and bottom-quark mass determination. This method 
has first been suggested in Ref. [1, 2] and has been continuously improved over the 
years by considering new data from experiment and by including higher order correc- 
tions [3-5] . However, also the scalar, pseudo-scalar and axial- vector current correlators 
can be employed in this context. For example, moments of the pseudo-scalar correlator 
can be used in combination with lattice calculations to determine fundamental con- 
stants of QCD, the strong coupling constant and the charm-quark mass with high 
accuracy [6]. Furthermore, the first eight three-loop moments have been used [7-10] as 
one of several ingredients in order to reconstruct the complete momentum and mass 
dependence of the scalar, pseudo-scalar, vector and axial-vector polarization functions 
by means of Pade- approximations. The complete mass and momentum dependence is 
important in the computation of the Z-boson decays, in which a combination of the 
vector and axial- vector density enters or in the calculation of Higgs-boson decays which 
are related to the scalar or pseudo-scalar correlator. 

Recently the low energy expansion at three-loop order has been extended and the 
first 30 coefficients have been determined for the vector current and for the remain- 
ing scalar, pseudo-scalar and axial-vector currents in Ref. [11,12], where also singlet 
contributions have been taken into account. In four- loop order only the first moment 
of the vector current correlator is fully known at present [13, 14]. For double fermionic 
contributions the first five moments have been computed [15]. Contributions of the 
order ain{~^ have been calculated to all orders j in Ref. [16], where ni denotes the 
number of light quarks, considered as massless. Numerical results for two moments of 
the pseudo-scalar correlator were presented in [6]. 

Whereas the lowest moments of the vector correlator have been studied up to four- 
loop order already some time ago, the ones of the scalar, pseudo-scalar and axial- vector 
current correlators were available up to three-loop order only. The purpose of this paper 
is to provide the still unknown four-loop QCD corrections to the lowest moments of the 
scalar and axial- vector current correlators and to present the details of the evaluation 
for the pseudo-scalar correlator, where first numerical results were presented in [6]. 
Our discussion will be hmited to the non-singlet contributions. These results can be 
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useful in combination with lattice calculations to reduce the error due to unknown 
higher order corrections in perturbation theory in the context of the determination of 
the strong coupling constant and quark mass as demonstrated in [6] . They can also be 
seen as a first step towards the evaluation of higher moments. Furthermore analytical 
results of the pseudo-scalar correlator are presented. 

The techniques used in this work have already been successfully applied in sev- 
eral other calculations, among which are the calculation of the matching relation of 
the strong coupling constant at a heavy quark threshold up to four-loop order in per- 
turbative QCD [17, 18] and the computation of the four-loop QCD corrections to the 
p-parameter arising from top- and bottom-quark loops [19-22]. 

The outline of this paper is as follows: In Section [2] we introduce our notations and 
conventions. In Section [3] we discuss the methods of calculation and give the results for 
the first moment of the scalar and axial-vector current correlators at four-loop order. 
For completeness we recall the results for the pseudo-scalar and vector case. Finally 
in Section H] we close with a brief summary and our conclusions. The results for the 
vector and pseudo-scalar correlator are listed in Appendix \^ those for the moments 
with n = —1 and n = in Appendix [B] 

2 Generalities and Notation 

The polarization functions for the scalar (s), pseudo-scalar (p), axial- vector (a) and vector (t>) 
current correlator are defined by 

q^U\q^) = iJdxe''^''{0\Tj\x)j\0)\0), for5 = s,p (1) 
iq,q. - q'9,u) n^(g') + g^^.n^g^) = ^ Jdxe^^^ {0\Tf^ix)fM\0) , for 5 = a,v (2) 
with the currents 

The low-energy expansion of the polarization functions in 2; = g^/ (2 m)^ is conveniently 
written as 

nV) = T7^ E c^n^^ (3) 

where the expansion coefficients are computed up to four-loop order in perturbative 
QCD. The expansion in the coupling constant Os/tt up to four-loop order is given by 

7^ = + (^) ci."'' + {^yc^:'' + {^y + . . . . (4) 
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We decompose the coefficient ' (i = 0, 1, 2, 3) for n > into the non- logarithmic 
and logarithmic parts 

c^' = tc^'''iL (5) 

j=0 

with Im = log (^)- Furthermore we classify the diagrams with respect to the number 
of closed fermion loops inserted into a diagram. The symbol ni denotes the number 
of light quarks and the symbol Uh = 1 denotes a heavy quark with mass m. This 
decomposition at four-loop order is given by 

C^J^' = c[^^' + Hh + c[^f' rii + cl^^hh ^\ + ^hni + C^j/' n]. (6) 

The bar indicates that renormalization of m, as and the current has been performed 
in the MS-scheme. We have checked that for n = 1 both scalar and axial-vector po- 
larization functions 11 (g^) obey the standard renormalization group equation(RGE). 
For the vector current correlator the longitudinal part Tl\{q^) of the polarization func- 
tion is zero due to the vector Ward-identity. The longitudinal part of the axial-vector 
correlator obeys the axial Ward-identity [23-27] 

g4 Ul{q^) = (2 m)2 H'^iq^) + contact term. (7) 

Inserting the expansion of Eq. ([3]) leads to 

oo oo 1 

E Cln = E C^n + - contact term. (8) 

n=-l n=-l 1 

Performing a shift in the summation index n = k + 1 on the r.h.s. of Eq.(l8]) and 
comparing the coefficients of the different orders in in both sides leads to 

CL,k = C'fe+i, (9) 

for k > —1, which allows to obtain the second moment of the pseudo-scalar correlator 
from the calculation of the ffist moment of the longitudinal part of the axial-vector 
correlator. The contact term in Eq.([8]) only contributes to the order l/z"^. 

3 Calculations and Results 

In a ffist step the program QGRAF [28] has been used to generate all necessary diagrams. 
Subsequently all appearing integrals have been mapped on a small set of 13 master 
integrals with the traditional Integration-By-Parts(IBP) method [29] in combination 
with Laporta's algorithm [30,31]. This procedure has been coded with the help of the 
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programs FORM [32-34] and FERMAT [35]. The remaining master integrals are shown in 
Fig. [1] in the standard basis. They have first been determined in Ref. [36] with the 
method of difference equation [31,37] and subsequently in Ref. [38], where the method 
of e-finite basis has been introduced. Also other authors [19, 39-44] have contributed 
in this connection with analytical or numerical results. 




-■^71 ^72 -"-81 -"-gi 

Figure 1: Master integrals in the standard basis. The solid (dashed) lines denote massive 
(massless) propagators. 



The four-loop moment with n = 1 is known analytically for all four correlators. For 
the scalar correlator it is given by 
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2009 4 1937539 , 

H TT^H , (14) 

186624 4199040' ^ ' 

Ms _ 7867 

~ 32805' ^ ^ 

W3i),s _ 197329 2573 / 2541989 _ 429089 

^ ~ 1296 ~ ~W ^ V 233280 ~ 31104 
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/ 1597 1435 \ 2 238 
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Ci = Uh rii \- rih hninh h n, , (17) 

^ 1620 1215 1215 '^3645 3645 '3645' ^ ^ 

where denotes the Riemann zeta- function and a„ = Li„(l/2). The result for the 
axial- vector correlator is given by 

-=(30) a 3996704 5966100779 499588, 
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59917211 ^ \ /22748503 14152979 , , 
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Numerical results for the pseudo-scalar case have first been presented in Ref. [6]. The 
analytical result is given in Appendix |Al together with the corresponding one for the 
vector current correlator, taken from Ref. [13,14]. The coefficients {5 = s,a,p,v) 
and are listed in numerical form in Table [H 



s 


n 


■pW.d -^(10),d -^(20),A -^(21),A -^(22,d) -^(30),A -^(31),A -^(32),d -^(33),A 


s 


1 


0.8000 0.6025 0.0000 -7.7402 -1.2551 0.0000 -5.4135 32.8284 2.6149 0.0000 


a 


1 


0.5333 0.8461 1.0667 -1.1071 1.6925 -0.0444 -2.4297 5.2981 0.4406 0.0321 


P 
P 


1 
2 


1.3333 3.1111 0.0000 0.1154 -6.4815 0.0000 -1.2224 2.5008 13.5031 0.0000 
0.5333 2.0642 1.0667 7.2362 1.5909 -0.0444 7.0659 -7.5852 0.5505 0.0321 


V 


1 


1.0667 2.5547 2.1333 2.4967 3.3130 -0.0889 -5.6404 4.0669 0.9590 0.0642 



Table 1: Numerical values for {6 = s,a,p,v) and for n/ = 3 which corresponds 
to the case of the charmed quark. 

In Appendix |B] we provide in addition the expansion coefficients for n = 0, —1 . Except 
for Cq they are only known numerically, however, with high precision. A completely 
analytical result would require the analytical determination of the constants T^4,2, 764,2, 
^61,2, ^"62,3, 772,1, T?!,!, T'gi,!, Tgi i, where Tn^i are the coefficients of the e-expansion 
(e = (4 — d)/2) of the master integrals T„ as defined in Fig. [T] 

oo 

T„= ^ e'Tn,i. (27) 

4 Summary and Conclusion 

We have computed the lowest coefficients in the low-energy expansion of the scalar and 
axial- vector current correlators at four-loop order in perturbative QCD. All appearing 
loop-integrals have been reduced to known master integrals, using Laporta's algorithm. 
We also gave the details of the calculation as well as the analytical results for the mo- 
ments of the pseudo-scalar correlator. These results allow to reduce the theoretical 
error originating from higher order corrections in the determination of fundamental 
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constants of QCD, like the strong coupling constant and the charm-quark mass in the 
context of lattice calculations. 
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A Moments for the pseudo-scalar and vector 
correlator 

The analytical result for the first two moments of the pseudo-scalar correlator presented 
in numerical from in Ref. [6] is given by 
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The first moment of the vector current correlator has been obtained in Ref. [13, 14] 
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348701 , 2/ N 2 1447057 . 27670774337 
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B Moments for n = — 1, 

The moments for n = —1 and n = exhibit an overall UV divergence. The correspond- 
ing l/£-poles are dropped by definition in the MS-scheme. The finite parts are given 
by 

= -325.6276432 + 16.39537650 Uh + 19.76434509 m 

-1.670198265 - 0.9856898698 n^rift + 0.7103788267 nf, (54) 
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+0.08246295965 nl - 0.5283433562 ni uh - 0.4972293123 , (57) 
C^ll" = 25.62696915 - 8.147150256n^- 6.938402913 
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-0.08246295965 + 0.5283433562 rii rih + 0.4972293123 nf, (58) 
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where for compactness /i — m has been chosen. The logarithms log(/x^/m^) can be 
reconstructed with the help of the RGE. 
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